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Abstract 

A simple model of macroevolution is proposed exhibiting both the 
property of punctuated equilibrium and the dynamics of potentialities for 
different species to evolve towards increasingly higher complexity. It is 
based on the phenomenon of fractal evolution which has been shown to 
constitute a fundamental property of nonlinear discretized systems with 
one memory- or random-based feedback loop. The latter involves a ba- 
sic "cognitive" function of each species given by the power of distinc- 
tion of states within some predefined resolution. The introduction of a 
realistic background noise limiting the range of the feedback operation 
yields a pattern signature in fitness space with a distribution of tempo- 
ral boost/mutation distances according to a randomized devil's staircase 
function. 

Introducing a further level in the hierarchy of the system's rules, the 
possibility of an adaptive evolutionary change of the resolution itself is 
implemented, thereby providing a time-dependent measure of the species' 
cognitive abilities: an additional feedback loop makes use of the inevitable 
intrinsic fiuctuations within the fitness landscape to direct the temporal 
change of the resolution. Feeding back the small adaptive changes of reso- 
lution into the essentially directionless variations of the patterns' lifetimes 
in fitness space effectively leads to a clear tendency towards increasing 
evolution potentials for each species ("hierarchically emergent fractal evo- 
lution"). 

Keywords: Macroevolution, punctuated equilibrium, fractal evolution, 
growth of complexity 
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1 Introduction: Punctuated equilibrium and self- 
-organized criticality 



Punctuated equilibrium is a well-documented extension of Darwinian evolution- 
ary theory On a population level, it has recently been observed in 
the laboratory: during its evolution, the average cell size of Escherichia coli 
bacteria showed long periods of stasis with rare beneficial mutations rapidly 
sweeping through the whole population Presently, however, it is not 
yet clear, whether punctuations in the domain of populations and species, re- 
spectively, are unrelated, or consequences of the same evolutionary mechanisms 
at work. On the species level, some of the qualitative features of punctuated 
equilibrium have been successfully simulated by Bak and Sneppen in a toy 
model on a lattice of sites, with each site representing a species, thereby drawing 
on earlier work on self-organized criticality (SOC) Q. Recently, exact results 
have been obtained for spatiotemporal evolutions in an extended SOC model of 
punctuated equilibrium with many internal degrees of freedom per site Q , for a 
model with neural- like rules ||l9| , and for generating phylogenetic-like trees |^ . 
However, the successful implementation of SOC models notwithstanding, there 
exists one important question which is hardly touched upon by these models, 
namely how evolution produces entities of increased complexity. In particular, 
it has been proposed by Kauffman and coworkers that the Darwinian process 
of mutation and environmental selection would have to be complemented by 
self-organizing processes in order to eventually decouple the strict link between 
mutation and selection p^ . 

In this paper we present a simple toy-model of evolution which is general 
enough to be compatible with SOC models of punctuated equilibrium, but in es- 
sential ways goes beyond them to encompass evolution towards higher complex- 
ity. For, if one considers the Bak-Sneppen model and its more recent versions, 
one is led to conclude that the SOC state is always a "quasi - steady state", i.e., 
"the evolution process always self-organizes into the same critical steady state 
having, as a consequence, the same appearance" [Q. In other words, the anal- 
ogy to the well-known sandpile model suggests that SOC describes evolution via 
punctuated avalanches of high mutation rates over a wide range of sizes, with 
long intermittent periods of stability, but without any (long term) time depen- 
dence of the avalanche sizes, the equilibrium periods, or their relative changes, 
respectively. This overall large-scale time symmetry is inherent in the model of 
Bak and Sneppen, as is easily seen from its results. 

In their one-dimensional case that interests us here, each site on a lattice 
representing a species interacts with its two nearest neighbors (thus representing 
a food chain, for example). The individual species (sites) are assigned random 
numbers between and 1 representing their fitness in the fitness landscape of the 
whole evolving pattern of lattice sites. At each time step, the site with the lowest 
fitness is chosen to be replaced by a different random number, thus mimicking 
Darwinian extinction of the least fit species and its replacement with another 
species in the same ecological niche. The same procedure is also applied to the 
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two nearest neighbors of the least fit species to simulate coevolution, so that in 
sum the fitness values of three sites are changed per time step. The result of this 
procedure is that most species self-organize into uniform distribution above some 
critical fitness value, while the few below it produce said punctuated avalanches 
of various sizes. The resulting activity pattern for the mutations clearly shows 
punctuated equilibrium behavior which is time-symmetric over long periods of 
evolutionary time (cf. Fig. 6 in |^). Moreover, as the model allows to add 
together only accumulated mutations at single sites, i.e., without differentiating 
between progress and regression in fitness over time, the accumulated activity 
pattern is also time-invariant (cf. Fig. 7 in ||] ) , and has recently been identified 
with the fractal pattern of a randomized devil's staircase function 

In contrast, we design our toy model in such a way that four major problems 
are directly adressed to overcome some weak points in the existing models of 
punctuated equilibrium: 1) To potentially allow for the emergence of ever higher 
complexity, the implicit overall time-symmetry of the evolving patterns in state- 
space must be discarded. 2) Accordingly, a plot of evolutionary changes for a 
single species over time must not only add accumulated changes (producing a 
devil's staircase function), but must be able to differentiate between progress 
and regression with regard to some measure of fitness or complexity (cf. Fig. 

1 in and Fig. lA in 3) The strict coupling between mutation and 
selection, i.e. the ongoing and strictly time-ordered interplay of both mecha- 
nisms, which is often considered as a necessary prerequisite for evolution, has 
to be relaxed to allow for the self-organization of evolutionary patterns. In this 
way, coevolutionary processes can lead to an effective decoupling of mutation 
and selection. 4) As the model to be designed should depend even less on the 
concrete content of micro-evolution, it should be as abstract as possible to study 
potential trends in evolution with a high chance for universality, i.e., irrespective 
of the underlying detailed processes. We thus aim at a meta-macroevolutionary 
model in the sense that we disregard questions on the concrete evolutionary 
substrate, thereby concentrating on relative changes with respect to previous 
evolutions. In other words, we are interested in simulating the "evolution of 
evolution". Also, as far as the assumption holds that the same evolutionary 
mechanisms are at work at the levels of populations and species, respectively, 
our model for the species level will apply to the population level as well. 

To implement the four points just mentioned, we propose the following simple 
toy-model for the evolution of evolution. 

2 Hierarchically emergent fractal evolution (HEFE) 
with noise 

2.1 Fractal evolution 

We consider the evolution in time of a one-dimensional array of cells poit,j) 
with j — I, . . . , N representing the number of species on a coupled map lattice 
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(CML). The nearest neighbor interaction implemented is of the form 



Poi*^^) = mtjiPoi* - 1. j) + -5 • [Po{t - 1, J - 1) + po{t - 1, J + 1)]} 
j^l,...,N, po,5en 

with S as couphng constant and A/q as normahzation factor yielding 

N 

J2po{t,j) = l foraUi. (2) 

i=i 

To constrain the values of the sites within a realistic domain, we introduce 
a system's internal precision hmit ^max = 10^® with the effect that each site 
Poit^j) < £,max is replaced by a white noise term, generated by ^max ■ RAN. 
The random number RAN is uniformly distributed on the interval [0, 1]. 

In this way, po{t) assumes values between zero and one representing the 
individual species' fitness at time step t relative to the others in such a manner 
that the total fitness for the whole array is given by a constant which, without 
loss of generality, is defined to be equal to unity. In other words, just as in we 
assume local interaction between the nearest neighbored species to affect relative 
fitness. However, as opposed to we are not interested in implementing a 
very specific and thus extremely model-dependent mechanism of eliminating 
least fit species, or the like. Rather, we intend to study the relative frequency of 
deviations from some average evolution, the details of which are not supposed 
to matter. To do so, we need to have an evolution of the form of Eq. (1) upon 
which certain constraints can be imposed, such that some kind of statistically 
controllable but individually contingent deviation can be implemented. 

With regard to this, a phenomenon is at hand which has been discovered and 
studied by us in some detail, called fractal evolution ^, . It is a universal dy- 
namic property shown for one-dimensional CMLs with (i) one or more temporal 
feedback operations (involving some memory of the system's states) boosting 
certain sites to values far above the average ones, and (ii) a normalization pro- 
cedure after each time step as given by Eq. (2). Fractal evolution is character- 
ized by a power-law behavior of the system's order parameter with regard to a 
resolution-like parameter which controls the deviation from an undisturbed (i.e., 
feedback-less) system's evolution and provides a dynamically invariant measure 
for the emerging spatiotemporal patterns [||. It must be strictly distinguished 
from the phenomenon of fractal growth, which describes the accumulation of 
micro-patterns into a static structure with (self-similar or other scale invariant) 
fractal features. In contrast, fractal evolution is characterized by the fact that 
the dynamics of the evolution itself is scale invariant: instead of generating a 
"frozen" fractal object, the generation mechanism itself is "fractal", i.e. the cho- 
sen resolution itself finally generates the observed fractal properties. In reference 
[ p^ we showed that the phenomenon is very robust with respect to variation of 
the system's variables and holds also if, instead of a temporal feedback condi- 
tion, the deviations from the average evolution are implemented by randomly 
generating local boosts to high values for individual sites. In these cases, the 
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only difference to the feedback situation is a different resulting fractal evolu- 
tion exponent in the power law. This means that no matter how the detailed 
micro-evolutionary substrate is implemented, i.e. by either entirely determinis- 
tic temporal feedback loops with a memory, or by completely random boosts to 
produce patterns deviating from the average development of the system, fractal 
evolution will be observed in most cases pO[ |. 

In the following, we shall choose temporal feedback loops in our model be- 
cause they introduce a quality of self-reference for the whole system which can 
be expanded upon later (cf. Sec. 2.2). As a concrete motivation for doing 
so, we consider a certain periodicity of constraining conditions for the whole 
evolutionary landscape. Again, the details are not supposed to matter, i.e., 
the periodicity may refer to relatively short-term cycles of constraints enforcing 
functional couplings within individual species, or externally, in the ecology or 
the like. However, the period of constraints should also not be too short, be- 
cause otherwise internal fluctuations would merely produce spurious effects of 
the internal systems dynamics which are of no particular interest here. 

In a "normal" evolution of our landscape that were unaffected by such cy- 
cles, the relative fitness values would gradually become adjusted to each other, 
thus gradually reaching an ever more smeared-out pattern representing stasis or 
equilibrium. However, if a species' fitness value after such a cycle happens to be 
the same again as before it, within a certain degree of accuracy 1/e, then, even 
if the relative fitness value is low, some reproductive isolation must have been 
achieved (for example, via economic functional couplings) that amounts to an 
acquisition of a relative survival rate far above the average. In fact, "molecu- 
lar cognition" within genomic dynamics [ pT[ already suggests a certain cognitive 
ability for each individual site in that it can monitor its own "fitness" within that 
of its environment with a resolution 1/e. Maintaining constancy in fitness space 
within cycle periods with a constraining and eventually destructive potential, 
while most of the environment produces decreasing fitness values, means that 
an ecological niche has been found whose fitness potential can be consumed. 
(The overall relative fitness is thus decoupled from the individual and momen- 
tary survival rate: the former is a purely coevolutionary measure, while the 
latter refers to a concrete situation in the larger context of cyclically changing 
constraining conditions of the whole evolutionary landscape.) In these cases, 
then, the species' site is a particularly favorable niche and therefore enhanced 
to some boost value, thus providing a fragmented overall pattern emerging in 
fitness space with particular elevated zones of highest survival rates (and, con- 
sequentially, also of fitness). 

These zones represent a high possibility for reproductive isolation (for ex- 
ample on the population level) and thus also for protected change. Although in 
principle it can occur at any time or site, a species' possible mutation, then, is 
most likely to accumulate into an evolutionarily conserved change (for example 
a morphological one) if occuring in such boosted niches of survival rates (and 
thus also new fitness values) far above the average. Naturally, if the accuracy 
1/e within which present and past relative fitness values are observed (i.e., to 
produce constancy with respect to the steadily evolving environment) is rather 
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low, there will be a relatively high chance for various boosts to occur within 
short periods of time. If, on the other hand, 1/e becomes higher and higher, 
i.e., if the cognitive requirements regarding the resolution become ever more 
demanding, then the chance for boosts gradually diminishes, and the actually 
boosted patterns will survive for longer periods of time. 

To implement said temporal feedback operation in our model of Eq. (1), we 
install an enhancement (of strength Aamp) of each site p{t,j) > ^max fulfilling 
an e-condition 



with 

Po,dtJ) = [Po(t, j) - • % + C , PoitJ) defined in Eq. (1) , 

e = Uax ■ RAN, RAN e [0, 1], = e{po{t,j) - imax) , 



e, = e 



P(t-t™em,i) 

and @{x 



<€<!{= 100%) , emem = &{t - tmem) 

1 : a; > 
: X < . 



(3) 

As usual, we impose periodic boundary conditions. Also, we introduce a time 
span tmem during which the system's memory is effective. (For example, the 
latter may reflect the maintainance of a particular molecular state during one 
cycle of the periodical constraining conditions.) 

In other words, the enhancement for any site Po,(^{t,j) > £,max takes place if 
t > tmem (i-c, the time span during which a system's memory is effective) and 
if the condition p{t - tmem J) • (1 - e) < po.d^J) < pit - tmem J) • (1 + e) is 
fulfilled. Note that this kind of feedback operation simulates a selection process 
enhancing the values and, consequently, the further evolutions of those sites 
which are comparable in fitness with their corresponding values in the past 
since they fulfill the e-condition. 

The whole normalized system is constrained to the range < p{t,j) < 1. 
Values of p{t, j) < ^max are, just for simplicity, replaced by white noise terms, 
i.e., in that case, whatever kind of noise existing in our system exceeds a signal 
generated by the system's rules. Moreover, a further parameter is introduced 
as an external representation threshold L » ^max- with it a background is 
defined against which the arising patterns obeying p{t,j) > L can easily be 
discerned. Otherwise, e.g. for L = ^max, the evolving systems' pattern would 
appear as one connected stream and could not be discerned into the separated 
fragments of the landscape's peaks needed for quantitative analysis. 

For the further analysis and without loss of generality (cf. the discussion 
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Figure 1: Temporal evolution of an arbitrary site (tt48), initiating from the noise 
domain below the threshold £,max, and later entering the region of an observable 
pattern above the representation threshold L. 



below) we choose a default set of parameter values: 
Number of sites (dimension) = 120 

Initial values (before normahzation) p{0, 40) = 0.1 , p(0, 60) = 0.9 , p(0, 100) = 
Mixing parameter (coupling constant) 5 = 0.03 
Boost value (enhancement) Aamp = 100 
Size of memory tmem — 200 
Threshold value L = 10"^ 

(4) 

In Fig. (P the time evolution due to Eq. ^ of one arbitrarily chosen site 
(i.e., site (148 in our case) is displayed together with the noise threshold (dotted 
line) and the external representation threshold used in the further graphical 
representation of the plots (dashed line). Due to Eq. (||), at time steps t — 1046 
and t = 1108 the site is boosted up to a value of order 0(1) after normalization, 
whereupon the typical decrease due to the diffusion rule (|l|) is observed until 
the boost of some other site within the array suppresses it by the factor Aa,np = 
100 due to the normalization procedure. The feature of repeated staircase-like 



7 



Figure 2: Density plot of the evolving system accentuating the range of the box 
values from L = 10~^ to 10""'^. Values larger than 10~^ are displayed with the 
maximal box size. The relative interval width e is chosen as e = 30%. 

increases of the site's values comes from contributions of neighbored large- valued 
sites. 

A spatiotemporal plot of all sites above the threshold value L is shown in 
Fig. @). Contrary to the signature of fractal evolution introduced in |l^, the 
fragments gather around a narrow array and are not distributed over the whole 
spatiotemporal plane. This effect represents a certain "taming" of evolutionary 
chance and comes from the restriction requiring that all sites be larger than the 
noise threshold £,max for their possible enhancement. (As stated above, eventual 
values below the noise threshold thus contribute only to the noise itself, but 
cannot be boosted to values dominating the whole evolution.) Since the range 
of those values which are checked by the e-condition of (||) contains only six 
orders of magnitude, a boosted site will be placed in the spatial neighborhood 
of a fragment tmem time steps ago. As will be shown below, this "stream- 
like" feature in the fitness landscape is crucial for the punctuated equilibrium 
behavior during the systems evolution. 
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Figure 3: Log- log plot of the mean lifetimes f of the fragments versus the 
relative interval width e. System parameters are taken from Eq. (^). The 
fractal evolution exponent is estimated by the slope of the fit curve. 
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Figure 4: Same as in Fig. but now only for each 1000th time step, in 

order to see the long-term evolution. The narrow band of fragments in Fig. 

eventually enters practically all regions of the array in a random-walk-like 
manner. 

In Fig. (H) the system's global property of fractal evolution is shown. The 
mean lifetime f of the patterns which is estimated by the arithmetic mean of 
the maximal temporal extension of all fragments for a chosen e (cf. Fig. (||)) 
scales according to 

f = a•e^ with a = 2.69 ±0.05, 5 =-0.620 ±0.005 , (5) 

practically irrespective of variations of the systems variables or the initial con- 
ditions. The power-law behavior of the system's order parameter f with regard 
to the resolution-like parameter e provides a dynamically invariant measure for 
the emerging spatiotemporal patterns given by the fractal evolution exponent 
b. As will be discussed in Sec. 2.2, fractal evolution in its dynamic sense is 
obtained only if the restriction to one specific value of e is abolished. Note that 
the system evolving under a chosen e does not evolve automatically towards 
a critical state in the sense of the SOC theory Q since the avalanches of the 
fragments' sizes do not occur in our system on all scales. 
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Although the two concurring conditions in Eq. (Q), i.e. the temporal e- 
condition and the noise-threshold condition, constrain the evolving patterns 
to lie within a narrow range, the long-term behavior as shown in Fig. 
exhibits a random-walk-like movement of the stream of fragments across the 
whole array. The invariance of this ergodic-like property under variation of 
the systems parameters has been checked empirically. Even the replacement 
of white noise values for sites below £^raax by a constant value (by 10~^, for 
example), provides the observed ergodicity, although the stream of fragments 
then evolves symmetrically due to the totally symmetric evolution rule and the 
uniform background. Therefore, in the long run, practically each site/species 
gets the chance to obtain a high survival rate and, consequently, to evolve in a 
high fitness potential. 

Any enhancement of a site represents a strong intervention for the system 
as a whole (due to the normalization) and for the represented species itself, 
of course, since it can be interpreted both as a drastic change in coevolution- 
ary fitness or as a site-bound beneficial mutation, respectively. Consequently, 
the time span between two succeeding boosts for any site offers both the op- 
portunity for further beneficial selection during its evolution, or towards ever 
decreasing fitness, and ultimately being swept below the noise limit ("detrimen- 
tal extinction"). In Fig. (H) the accumulated boosts are displayed for a time 
interval of 50000 time steps. Similar to the SOC models |§, |^, each site exhibits 
a punctuated equilibrium behavior with respect to the time plateaus indicat- 
ing the time difference between two drastic evolutionary activities. Due to the 
chosen intial conditions and within the first 50000 time steps, site pO happens 
to experience many more changes than site JJIO. Since the stream of fragments 
and thus the main boost activities will reach also previously quiescent regions 
after some time (cf. Fig. (^), the accumulated boosts of site jjlO will even- 
tually catch up with those of site pO. The dashed curve in Fig. (||) shows 
the effect of omitting only one boost (i.e., the 100th one in our case) during 
the site's evolution. The observed dramatic change of the accumulation curve 
demonstrates just one example of the crucial role of contingency in the systems 
dynamics. The analysis of the distributions of time plateaus for an arbitrarily 
chosen site yields a (discretized) shape as shown in Fig. (^. Each bin contains 
the relative frequency of ten succeeding plateau sizes. Apparently, practically 
all sizes of periods of stasis are present, the data for bin p. up to about pO 
even lying on a staight line in a log-log plot. The peaks around bin pO, i.e. for 
plateau sizes of order 0(200), and another, weaker one around bin pO (plateau 
sizes of order 0(400)) represent a direct effect of the chosen size of the system's 
memory tmem = 200. Succeeding boosts occurring with temporal distances of 
order 0{k ■ tmem)i k = 1,2,... are causally related due to the evolution rule 
and the values of the systems variables. The remaining boost distances are 
clearly not directly related and their distribution represents a genuine emer- 
gent property of the systems dynamics. As long as we use a systems memory, 
this perturbation of causally unrelated boost distances by causally related ones 
seems to be inevitable, but the moment of the peak's first occurrence can be 
shifted towards larger plateau sizes as shown in Fig. m). In that case, a memory 
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Figure 5: The accumulated boosts for two different sites show punctuated equi- 
librimn behavior, i.e., periods of stasis arc interrupted by high activity of boosts 
including eventual mutations. The dashed line shows the effect of omitting the 
100th boost for site pO. 
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Figure 6: Relative frequency distribution of the sizes of the time plateaus in 
Fig. (^) according to their size, represented by bin numbers with a bin width of 
10. Bin jJlO, e.g., contains all sizes from 90 to 99. Bin jJlOOO contains the sizes 
9990 to 9999 and, additionally, all other ones larger than 9999. The data were 
taken over a period of 2.10^ time steps. Note the appearance of the first peak 
around bin (J20 which is explained in the text. The left hand side of the figure 
represents behavior according to a randomized devil's staircase function. 



13 



Figure 7: Same as in Fig. (^), but for a larger memory tmem = 800 and for 
a larger number of bins to make the occurrance of very large time plateaus 
more clear. The first peak now shows up at around bin [180. The circled data 
represent the originally introduced noise-free model of fractal evolution ||, 
indicating a sharp cutoff of the plateau sizes. The fit indicates the slope due to 
a randomized devil's staircase function. 
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depth of tmem = 800 was chosen. Accordingly, the peak shows up at about bin 
U80. Up to this value, the frequency /(s) of plateau sizes s scales according to 

/(.) oc ,-i-i7±o.oi (g) 

thus indicating a randomized devil's staircase function |^]. The exponent 
obtained above represents a robust property of the system in spite of its inherent 
contingency indicated by Fig. (||). For comparison, the circled data in Fig. (0) 
represent the distribution of boost distances due to our evolution rules used 
earlier where no noise threshold condition had been incorporated. The 

observed fragments' distribution all over the spatioteniporal plane in the noise- 
free model is reflected in the sharp cutoff of frequencies for plateau sizes beyond 
s « 10^. 



2.2 Self-organized evolution of fractal evolution 

As has already been remarked in the previous section, fractal evolution in the 
true sense of the word is effectively realized only if the interval width e is allowed 
to acquire different values during the systems evolution. According to a recently 
introduced idea ||ll|, to be further elaborated here, we impose an additional 
long-term control cycle leading to an adaptive change of e on the next higher 
level of the system's hierarchical nesting of feedback loops. 

We thus let the system be driven by fluctuations of fexp measured with the 
same e during two succeeding long time spans or generations {n,n + 1). The 
difference fexpin+1) —fexpin) is practically always nonzero due to the finite time 
interval of data taking. To implement an adaptive change of e or resolution 1/e, 
respectively, with regard to an ever changing ruggedness of the whole fitness 
landscape, this difference is fed back into the consecutive value of e via the 
difference quotient of the power-law (^) 



T e 



leading to 



. £n Tn+l Tn / , 

Ae = e„+i - e„ = . (8) 

Tn 



Since the last term of Eq. (g) is generally non-zero due to the fluctuations of 
the mean lifetime, a different relative interval width results governing the two 
succeding generations, and is determined by 

f. , 1 Tn+l -Tn\ , , 

e„+2 := En+l = En I 1 + ^ T 1 • (9) 

The geometrical analogy to the analytical determination of the next value of 
e„+2 or e„+i, respectively, is given as follows: the latter lies at the intersection 
of the tangent at (e„,f„) belonging to the theoretical curve f„ = a' • on one 
hand, and the horizontal line through fn+i on the other. 
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As a main result, we obtain a progressive alteration of e or the resolution 1 /e, 
respectively, within the system's evolution which we denote as "hierarchically 
emergent fractal evolution" (HEFE). In effect, HEFE constitutes a fluctuation 
driven and adaptive mechanism for the gradual increase of the mean lifetimes 
of the emerging patterns, i.e. for the expansions of regions of high fitness in the 
fitness landscape. 

In Fig. (^), three examples of temporal evolutions of the mean lifetime per 
generation are shown according to three different time spans of each generation 
during which the mean lifetime is determined. Although fluctuations are allowed 
with equal opportunity for decreasing and increasing lifetimes, respectively, a 
long-term increase of the mean lifetimes is observed practically in all cases and 
has its origin in the functional property of the power-law (^): Denoting the 
deviation of the mean lifetime by t„+i — f „ =: 77^ • f„, and assuming one 
positive-valued deviation 77"*" for the consecutive generations (n, n + 1) and a 
negative-valued one r]~ of about equal size w —77^ for the succeeding ones 
{n+2, n+3), one obtains altogether a clear tendency towards decreasing relative 
interval widths 



leading to an average increase of the mean lifetime f per generation. If f would 
depend, e.g., exponentially on e, like f cc e^'^, e would only oscillate between 
two fixed values. 

For relatively small time spans, i.e. for curve a) with 2000 time steps per 
generation, the fluctuation of the lifetimes is larger, of course, than for 10000 
time steps. Therefore, the observed drastic increase of the mean lifetimes will 
occur sooner. In sum, the systems dynamics makes possible a long-term growth 
of each species' lifetime and thus also of its complexity: larger lifetimes of high 
fitness periods inrease, e.g. via "learning", both the cognitive performances as 
well as the probabilty for survival of the species and, consequently, also the 
probability for genetic improvement via mutational mechanisms. As far as the 
cognitive functions are concerned, their improvement increases the lifetimes of 
high fitness values, and so on. We reach the limits of our model when the mean 
lifetimes in Fig. reach the size of the system's memory tmem , because then 
fractal evolution gets disturbed and new boundary conditions would have 
to be considered that might constrain the further evolution effectively. 

Finally, we briefly discuss the evolution of the system's entropy, which mea- 
sures the fitness of the whole ensemble of species, and is simply chosen as the 
relative Shannon information entropy |Ig] 



We can allocate the various activities within the systems evolution with 
regard to their specific entropy behavior in the following way. A site's enhance- 
ment, i.e., an increased single species' fitness value, lowers the system's entropy 
at that time step due to the enlargement of the values' differences within the 



(10) 




(11) 
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Figure 8: Hierarchically emergent fractal evolution (HEFE) for three different 
values of time spans of data taking with particular e ("generations"). The time 
span for each generation is given by 2000 time steps for curve a), 5000 for 
curve b), and 10000 for curve c). The starting value of e for all curves is 
eg = 10%. Apart from the obvious time asymmetry of the evolution, note also 
the characteristic stepwise increases as well as decreases of the mean lifetimes 
per generation, thus indicating two opposing trends in the short-term dynamics. 
In the long run, however, the mean lifetime is bound to increase. 
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Figure 9: Normalized entropy per 10000 time steps versus the relative interval 
width in a semi-logarithmic plot. Since the long-term behavior of HEFE devel- 
ops towards larger mean lifetimes f , and, therefore, towards smaller values of e, 
the entropy of the evolving system increases on the average. 

array. The effect of the noise threshold and the decay of the fragments due to 
Eq. (|l|) are entropy increasing processes. Since the systems evolution with a 
given fixed value of e yields a stationary pattern production characterized by 
the fractal evolution exponent, an average value of the above described entropy 
contributions summed up for an appropriate long time interval is obtained for 
that specific e. The larger the value of e, the shorter the mean lifetime of the 
fragments and, consequently, the more rugged the fitness landscape will ap- 
pear, thereby lowering the normalized entropy per time interval (cf. Fig. (^). 
With respect to the HEFE-mechanism, this means that the total entropy per 
time interval increases on the average due to the growth of the mean lifetimes. 
That is, within the scenario of macroevolution this effect can be considered as 
tendency towards maximal fitness of the whole ensemble of species under con- 
stant boundary conditions. Thus, internal evolutions within a species towards 
higher complexity are accompanied by a complementary, external increase of 



18 



the entropy of the fitness landscape. 



3 Conclusions 

The scahng behavior of certain observables in dynamical systems is of partic- 
ular relevance in a wide range of systems, be they geophysical, economical or 
biological, for example. Within the theoretical framework of the SOC theory 
some features of macroevolution have recently also become testable via a kind 
of "experimental paleontology" jl^. They are made plausible by the concept 
of a system's critical state which is reached after a certain evolutionary time, 
and which thereby exhibits a power-law behavior of the frequency of mutation 
activities with regard to their sizes. Yet they suffer from being trapped within 
the steady-state property of the system's dynamics. 

In this paper we have proposed a model for macroevolution based on concepts 
differing from the SOC model. Apart from some basic qualitative similarities 
(cf. Figs. and (||)), we obtain substantially different results, primarily with 
regard to the time-invariant behavior of the system's emerging order param- 
eter: the power-law obtained in our model reflects a universal feature named 
"fractal evolution", which can be based on a variety of underlying mechanisms, 
including the option of a memory-based feedback operation. The four major 
problems specified in the introduction in order to overcome some weak points in 
existing models of macroevolution are solved by the implementation of a second 
feedback operation. The latter transforms the mere functional property of "frac- 
tal evolution" of the system's fitness landscape via a fluctuation-based iterative 
mechanism into the dynamical process called "hierarchically emergent fractal 
evolution" (HEFE). In this way, a self-organized potential growth emerges on 
a long-term scale, such that the internal complexity of each species may grad- 
ually increase. However, throughout this work we have deliberately avoided a 
narrow definition of the term "complexity", or a choice among the proliferating 
definitions in the literature, respectively Instead, we have merely considered 
beneficial fitness domains for each species in the sense of potentials to evolve to- 
wards higher complexity, irrespective of the underlying concrete (self-organized) 
mechanisms. One main result of our study is that in our model such potentials 
grow with elapsing time. 

We also emphasize the important role of randomness in our model, which is 
effective on various hierarchical levels. At the basic level of pattern generation, 
we observe a strong sensitivity of the patterns' signature to slight variations 
of the initial conditions, similar to the features of deterministic chaos. An 
emerging global measure, irrespective of local system fluctuations, has been 
obtained in the form of the fractal evolution exponent. The incorporation of 
a noise threshold constraining the range of operational fitness values for each 
species enables the systems' patterns to walk randomly across the whole array 
(cf. Figs. (^ and (^), thereby producing activity /mutation patterns obeying 
a randomized devil's staircase function. As in the SOC model ^, small 
deviations in the order of mutation events produce significant deviations from 
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the course of the accumulated activity pattern (cf. Fig. (||)). The globally 
invariant measure in this case is given by the slope of the distribution function 
for the frequency of the periods of stasis versus their sizes (cf. Figs. (||) and 
(^). Moreover, at the level of the self-organized change of the resolution-like 
parameter e, the natural fluctuation of the order parameter f within two runs 
("generations") under the same conditions iteratively serves as driving force 
towards a changing e. The direction of this change appears on small time scales 
as a random walk. However, the long-time tendency towards the irreversible 
decrease of e, and thus towards the increase of the mean lifetimes of the high 
fitness domains for each species, respectively, is firmly anchored in the HEFE 
mechanism (cf. Fig. (||)). More generally speaking, it is rooted in the basic 
mathematical properties of the power-law describing fractal evolution, and it 
leads in the long run to a mutual positive enhancement of potential cognitive 
abilities for all species on one hand, and lifetimes with high fitness values on the 
other. 

Finally, the potential development of individual species towards gradually 
higher complexity (as represented by an evolution towards ever lower values of 
e) is accompanied by a tendency towards highest entropy values for the whole 
fitness landscape under constant boundary conditions (cf. Fig. (||)). The HEFE 
mechanism can thus be interpreted such that a single species is constantly en- 
gaged, via the possibilities of beneficial selection versus detrimental extinction, 
in a struggle for higher fitness values, whereas there exists a global tendency of 
ever higher fitness values to emerge for the whole ensemble of species. 
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